LU TP 96-26 
September 23rd, 1996 



MULTIPLICITY DISTRIBUTIONS AND FLUCTUATIONS IN 
THE DISCRETE QCD MODEL 

R. Ugoccionij] 

Department of Theoretical Physics, University of Lund, 
Solvegatan HA, S-22362 Lund, Sweden 



ABSTRACT 

Results arc presented from calculations carried out in the discrete 
QCD model, concerning multiplicity distributions in rapidity inter- 
vals and fluctuations. 



Talk presented at the 
7th International Workshop "Correlations and Fluctuation' 
(Nijmegcn, The Netherlands, June 30 - July 6, 1996) 



"E-mail: roberto@thep.lu.se 



MULTIPLICITY DISTRIBUTIONS AND FLUCTUATIONS IN 
THE DISCRETE QCD MODEL 



R. Ugoccioni 

Department of Theoretical Physics, University of Lund, 
Solvegatan 14A, S-22362 Lund, Sweden 

Results are presented from calculations carried out in the discrete QCD model, 
concerning multiplicity distributions in rapidity intervals and fluctuations. 



1 Introduction 



One of the long standing problems in multiparticle dynamics is an integrated 
description of multiplicity distributions and correlations. In particular, from 
the point of view of perturbative QCD, these observables are not satisfacto- 
rily described: this is particularly evident in restricted regions of phase space, 
where the most interesting experimental results areEl The discrete approxima- 
tion to QCDa allows to perform exclusive and inclusive calculations at parton 
level; in this talk, a few results are described concerning Multiplicity Distribu- 
tions (MD's) and fluctuations. 



2 A short description of discrete QCD 

Here I will present only the fundamental ideas behind this model; for a detailed 
discussion, follow the pointers in the bibliography!] Consider an event in which 
a quark-antiquark pair is produced: the cross section for the emission of gluons 
by this pair is, in the Double Log Approximation (DLA), given by the "colour 
dipole radiation" formulatl 

j 7} _ Ca s dk\ 

dP^—^dy (1) 

where C is a constant depending on the colour charges of the partons at the end 
points of the dipole, and y and k± are the rapidity and transverse momentum 
of the emitted gluon, respectively. Notice that the azimuthal direction is here 
integrated over, but it is of course possible to take it fully into account. The 
phase space limits for the gluon emission are given approximately by 

\y\ < log(W/fcx) (2) 

where W is the dipole mass (which in the case of the original qq dipole equals 
the cm. energy); the rapidity range available for the first gluon emission is 
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Figure 1: left) Phase space in the Dipole Cascade Model. The dimensions shown are y and 
k = log fc^ , respectively rapidity and log of transverse momentum, taken in each dipole's 
cm. frame, right) Phase space in the discrete QCD model. Emission of "effective gluons" 
can only take place at the points marked with a cross. 

thus Ay = \og(W 2 /k 2 J. After the first emission, the three partons form two 
independent dipoles: they will in turn emit gluons according to Eq. [I] (in 
the cm. frame of each dipole) and the cascade develops!] This is illustrated 
in Figure [l] (left part) where it is also shown that, because the gluons carry- 
colour charges, the phase space actually increases with each emission: a new 
double faced triangular region is added for each new gluon emitted; notice that 
according to Eq. || the baseline of each face added is half the height of the same 
face. The Modified Leading Log Approximation (MLLA) changes this picture 
in that it provides a region of phasei-space, close to the end points of the dipole, 
where gluon emission is suppressed!! Thus the effective rapidity range becomes 
smaller: for example, for a gg dipole one has for large W: 

W 2 

A y = lo § -p - ~ 5 V9 ( 3 ) 

It can be argued^'! from the form of the beta term in the Callan-Symanzik 
equation, that Sy g equals 11/6. Dipole splitting with ordering in k± can also 
be shown to bepjequivalent to the angular ordering prescription in standard 
parton cascades!™ _ 

In the discrete QCD modelp the idea expressed in Eq. is carried to its 
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extreme consequences: it is assumed that if two gluons are closer than Sy g they 
will effectively be re-absorbed into a single gluon. One needs then only consider 
"effective gluons" : the "rapidity" dimension of phase space is thus partitioned 
in discrete steps of size Sy g . However, this is true also for every triangle "added" 
to the original phase space, which means that discrete steps must be taken also 
in the vertical dimension, i.e., in k = log/c^: phase space is thus partitioned 
in rectangular cells of size 5y g x 2Sy g . In a complementary way, one can say 
that emission of effective gluons takes place only at a discrete lattice of points, 
Sy g units apart horizontally and 2Sy g units apart vertically (see Figure [j]). 
From the very structure of the cascade described above, it follows that each 
"column" of cells is independent of the others: such a column will be called a 
tree because of the possibility of emitting further triangles ( "branches" ) . The 
probability of emitting a gluon with transverse momentum k in such a cell 
can be obtained from Eq. [l] using the leading order expression for the running 
coupling constant a s and the appropriate colour factor for the g — > gg vertex: 



1 dn „ dn 
Sy g k 



dp = —— 5 y 9 = — ( 4 ) 



Normalizing with a Sudakov factor, the distribution in k for an effective gluon 
at a given rapidity is uniform: 

dProb = exp { - / dp\dP= (5) 



where K max corresponds to the maximum k± available at that rapidity. Each 
triangle emerging from the original one can then be described by the position 
and length of its baseline (both as integer multiples of Sy g ): an events is a 
collection of independent trees of different height, the number of trees being 
given by L — 2 \og(W/Sy g ). The maximum height of each tree (numbering, 
e.g., the trees from left to right) depends only on the tree's position j and on 
the event "baseline length" L. A tree of height H corresponds to an emitted 
gluon of transverse momentum k = H ■ 2Sy g . A tree of height corresponds 
to the emission of no gluons; notice that this introduces a natural cut-off for 
cascade evolution. Notice also that only a finite number of gluons can be 
emitted. These points are important in characterizing discrete QCD. 

In order to calculate the multiplicity distribution, one can now define an 
N -tree as a tree with height between and N — 1 (in units of Sy g ) where N 
is an integer. It is convenient to define also a true N -tree as a tree of height 
N — 1. Because of Eq. |, an iV-tree is obtained by summing true if-trees 
(H = 1, . . . , N) with equal probability 1/N. This is easily expressed in terms 
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of the generating functions for the MD's in an iV-tree, Vn(z), and in a true 
iV-tree, Q N {z): 



1 N 



H=l 



It will also be recognized that a true (N + l)-tree can be obtained from a true 
iV-tree by attaching two iV-trecs to it (one on each side) 

Qn+i(z)= Qn{z)[Vn{z)] 2 (7) 
Using the obvious initial conditions: 

Vt{z) = Qi(z) = l; Q 2 (z) = z; V 2 (z)= l -±^ (8) 

Eq.s ^ and (?] can be solved^ thus giving the MD in a tree. 



3 Calculations in rapidity intervals 

In order to approach the calculation in rapidity intervals, one should under- 
stand that the size of 5y g does not imply that there is no structure at distances 
smaller than 11/6. In fact, Sy g correspond to a separation in rapidity along 
the original qq axis (which one can identify to a good approximation with the 
thrust axis of the event) only for trees which 'grow' from the original baseline, 
while inside each branch, 5y g is a distance along the local 'jet' direction. Thus 
one can say that gluons belonging to the same tree are all localized in rapidity 
inside one unit of Sy g . That this is actually a good approximation has been 
checked with a Monte Carlo implementation of the discrete QCD model. In 
this framework one can then calculate observables integrated over intervals of 
rapidity whose width is a multiple of Sy g , i.e., which contain a whole number of 
trees. So one recognizes that the MD in any of these rapidity interval is given 
by V h l(z). Adjacent trees are independent, and so are adjacent rapidity in- 

3 

tervals in this approximation; thus the MD generating function in an interval 
Ay is given by 

V Ay (z) = [] Th-(z) (9) 

where the product is over all trees inside the interval Ay. 

Figure |L on the left, shows the resulting MD in different central rapidity 
intervals: the solid line shows a negative binomial distribution (NBD), phe- 
nomenologically widely used for describing the dataul The same figure shows 
on the right the MD for what I will call in the following "2-jet like" events, 
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Figure 2: a) Multiplicity distributions in central intervals of rapidity of size (from bottom 
to top) 2, 4, 6, 8 and 10 (in units of 5y g ) for standard events with baseline length L = 12. 
Each distribution is divided by 10 with respect to one above it in order to make the figure 
legible, b) Same as in a), but for "2-jet like" events. 



obtained by cutting away the tip of the phase space triangle, i.e., requiring a 
cut-off on the maximum k± of an emitted gluon. Notice that in the standard 
events, for small intervals, the NBD behaves w.r.t. the points in a way similar 
to the way it behaves w.r.t to experimental data on hadronic MD, going first 
below the data, than above and then below againQ In the 2-jet like events, on 
the other hand, the NBD is closer to the points than in the standard events: 
one is reminded of the experimental result that a NBD is found to describe the 
hadronic MD in. 2-jet events (defined via a jet-finder algorithm) better than in 
the full sample tl 



In the literature, the NBD is often explained in terms of clan structure 




independent emission of clans, according to a Poisson distribution, followed by 
particle emission from each clan in a cascade way (logarithmic distribution). 
In the framework of discrete QCD, it is tempting to notice that each tree which 
comes from the baseline of the original triangle is independent from the others; 
a tree which produces no partons is not "observed" , so the number of trees 
in a rapidity interval is variable. One can therefore call trees of height larger 
than non-empty trees; by definition these contain at least one parton. 
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Figure 3: a) Average number of non-empty trees in rapidity intervals vs size of the interval 
(in units of Sy g ), for three sizes of the event: L = 8 (crosses), L = 9 (diamonds) and L = 10 
(squares), b) Average multiplicity in an average non-empty tree vs size of the rapidity 
interval for the said event sizes, c) As in a), but for 2-jet like events, d) As in b), but for 

2-jet like events. 

Consider now an event with L trees: 

Xj = v H L(o) = -L (io) 

is the probability that the j-th tree is empty. The MD in the number of non- 
empty trees n is given by the product of the probabilities that L — n trees are 
empty and n are not: 

P n L) = J2^-xi)...(l-x n )x n+1 ...x L n = 0,...,L (11) 

perm 
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of the moments, q, for standard events of baseline length L - 

2-jet like events. 



vs the order 
10. b) Same as in a) but for 



where the sum is over all permutations of a;'s. It can be shown that Pn 
satisfies the recurrence relation 

^-E^!-^S (12) 

4=1 

where 




From these relations P n can be calculated; it is easy to show that one obtains 
a binomial distribution when Xj is independent of j, and a Poisson distribution 
when a r = 8\^ r : the former is incompatible with Eqj6|, but the latter is a good 
approximation when Xj is not very small. Figure |3| shows the result for the 
average number of non-empty trees N(Ay) and for the average multiplicity 
within an average non-empty tree n c (Ay) (defined as the ratio of the average 
multiplicity to N(Ay)) for standard events. The noticeable features for the 
average number of non-empty trees are an independence from the size L of the 
event in small intervals, and a linearity with the size of the interval at fixed 
L. Both these chaijaileristics are seen in experimental data with the statistical 
definition of clansLrEU 

Finally, in Figure || the results of a calculation in full phase space of the 
ratio H q of factorial cumulant moments to factorial moments are shown: these 
where shown to be sensitive to the structures in the MDjLHlO and in particular 
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to the radiation of hard gluons in the early stages of the partonic evolution. 
The figure shows that oscillations in sign of the ratio H q appear in the discrete 
QCD model, and when examined in the 2-jet like sample these oscillatkms 
diminish strongly in amplitude, in accordance with the behaviour of dataliij 
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